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fact,  the  results  show  anticyclonic  eddy  formation  during  periods  of  calm  Tehuantepec 
winds.  That  is  specifically  exemplified  by  the  formation  of  two  anticyclonic  Tehuantepec 
eddies  during  a  5-month  period  of  weak  Gulf  of  Tehuantepec  winds  during  summer  of 
1997.  Furthermore,  the  satellite-observed  and  NLOM-simulated  proliferation  of 
Tehuantepec  eddies  during  El  Nino  years  is  explained  by  the  corresponding  increase  in 
downwelling  coastally  trapped  waves  and  a  lack  of  increase  in  the  number  and  strength  of 
Tehuantepec  wind  events  during  El  Nino  years. 
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1.  Introduction 

[2]  Previous  studies  have  shown  that  the  generation  of 
oceanic  Tehuantepec  eddies  off  the  Mexican  southwest 
coast  is  due  to  intermittent  strong  offshore  winds  that  can 
exist  in  the  Gulf  of  Tehauntepec  (GT)  mainly  during  the 
boreal  cold  season  (fall -winter)  [ Clarke ,  1998;  McCreary 
et  al. ,  1989;  Lavin  et  al,  1992;  Giese  et  a /.,  1994;  Miiller- 
Karger  and  Fuentes-Yaco ,  2000],  The  GT  high  wind  events 
are  associated  with  the  arrival  of  high  pressure  following  the 
passage  of  atmospheric  cold  fronts  into  the  Gulf  of  Mexico 
f Barton  et  al.,  1993;  Trasvina  et  al.,  1995;  Bourassa  et  al . , 
1999;  Chelton  et  al.,  2000a].  This  high  atmospheric  pres¬ 
sure,  combined  with  low  atmospheric  pressure  in  the  Pacific 
Ocean,  creates  a  pressure  gradient  that  forces  southward 
winds  through  the  Isthmus  of  Tehuantepec  mountain  gap 
(Figure  1).  Those  winds  reach  the  Pacific  Ocean  at  the  GT 
with  gusts  of  35  m/s  [Romero-Cenleno  et  al,  2003],  As  a 
comparative  example,  a  tropical  storm  is  upgraded 
to  hurricane  status  when  its  sustained  winds  increase  to 
33  m/s. 

[}]  TOPEX/Poseidon  (T/P)  data  and  the  Naval  Research 
Laboratory  (NRL)  Layered  Ocean  Model  (NLOM)  results 
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show  interannual  variability  in  the  number  and  intensity  of 
the  Tehuantepec  eddies  (Figure  2).  Assuming  that  Tehuan¬ 
tepec  eddy  formation  and  modulation  are  predominantly 
due  to  the  GT  winds  we  would  expect  to  see  a  GT  wind 
related  interannual  variability  in  eddy  number  and 
strength.  However,  the  European  Centre  for  Medium- 
Range  Weather  Forecasts  ( ECMWF)  [1994]  (Figure  3) 
and  the  Fleet  Numerical  Meteorology  and  Oceanography 
Center  (FNMOC)  Navy  Operational  Global  Atmospheric 
Prediction  System  (NOGAPS)  [Rosmond  et  al,  2002]  (not 
shown)  wind  data  sets  do  not  show  the  hypothesized  GT 
wind  related  interannual  variability  in  the  Tehuantepec 
eddies.  In  contrast,  NLOM  simulations  forced  with 
ECMWF  and  NOGAPS  winds  show  that  besides  the  GT 
winds  (which  are  the  main  mechanism  for  Tehuantepec 
eddy  generation),  intraseasonal  (30-90  days)  and  interan¬ 
nual  (2-7  years)  equatorially  generated  poleward  traveling 
baroclinic  downwelling  coastally  trapped  waves  (CTWs) 
also  contribute  to  the  generation  and  modulation  of  the 
Tehuantepec  eddies.  In  similar  lines  of  research,  recent 
studies  have  discussed  and/or  invoked  the  role  of  CTWs  as 
eddy  generators  to  explain  eddy  formation  in  the  Gulfs  of 
Alaska  and  California  and  in  the  regions  close  to  Aca¬ 
pulco,  Cabo  Corrientes,  and  the  Maria  Islands  [Melsom  et 
al,  1999;  Murray  et  al ,  2001 ;  Zamudio  et  al,  2002,  2001 ; 
L.  Zamudio  et  al.,  Tropical  waves  induce  oceanic  eddies  at 
Cabo  Corrientes  and  the  Maria  Islands,  Mexico,  submitted 
to  Journal  of  Geophysical  Research ,  2004,  hereinafter 
referred  to  as  Zamudio  et  al.,  submitted  manuscript, 
2004].  Furthermore,  W.  S.  Kessler  and  Z.  Yu  (http:// 
www.pmel.noaa.gov/~kessler/t-peckers/phoenix-ams.html, 
1998)  using  an  ocean  general  circulation  model  noted  the 
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Figure  1.  TOPEX/Poseidon  (T/P)  sea  surface  height  (SSH)  anomaly  map  in  centimeters  for  the  period 
18-28  February  1993,  which  shows  the  presence  of  anticyclonic  Tehuantepec  eddies  and  their  associated 
surface  geostrophic  circulation  (arrow  vectors).  The  map  was  produced  by  interpolating  all  the  ascending 
and  descending  along-track  altimeter  measurements  to  a  1°  x  1°  regular  grid.  The  anomaly  is  with 
respect  to  the  1993-1999  mean.  The  positions  of  the  Gulf  of  Mexico  (GM),  Isthmus  of  Tehuantepec  (IT), 
and  Gulf  of  Tehuantepec  (GT)  are  indicated.  The  bold  black  line  identifies  the  satellite  track  along  which 
T/P  measured  the  Tehuantepec  eddies  shown  in  Figure  2. 


formation  of  Tehuantepec  eddies  when  GT  wind  events 
were  absent.  On  the  other  hand,  Hansen  and  Maul  [1991] 
proposed  an  alternative  generation  mechanism  for  the 
Tehuantepec  eddies  that  is  summarized  as  the  shedding 
of  eddies  by  the  North  Equatorial  Counter  Current  during 
its  change  of  direction  around  the  Costa  Rica  Dome. 

[4]  Here  we  investigate  the  origin  of  interannual  variabil¬ 
ity  in  Tehuantepec  eddy  generation  as  measured  by  satellite 
altimetry  and  simulated  by  a  numerical  model.  Two  hy¬ 
potheses  are  presented  and  discussed.  The  first  hypothesis  is 
that  the  GT  winds  are  the  only  generator  and  modulator  of 
the  Tehuantepec  eddies.  The  second  hypothesis  is  that  in 
addition  to  the  GT  winds,  intraseasonal  and  interannual 
CTWs  also  play  an  important  role  in  Tehuantepec  eddy 
formation  and  in  the  modulation  of  their  frequency  and 
amplitude.  Hypothesis  testing  includes  assessing  the  relative 
importance  of  GT  winds  and  CTWs  in  modulating  the 
interannual  variability'  of  Tehuantepec  eddy  generation. 

2.  Model 

[5]  The  primitive  equation  formulation  of  NLOM  has 
been  extensively  documented  by  Hurlburt  and  Thompson 
[1980],  Wallcraft  [1991],  Metzger  and  Hurlburt  [2001], 
Wallcrafi  et  al.  [2003],  and  Zamudio  et  al.  (submitted 
manuscript,  2004).  NLOM  application  to  Pacific  Ocean 
modeling  is  discussed  by  Hurlburt  et  al.  [1996].  The  main 
characteristics  of  the  core  NLOM  experiment  used  in  this 
study  are  summarized  below.  The  model  consists  of  7  layers 
(including  the  mixed  layer)  and  includes  a  free  surface, 


isopycnal  outcropping,  realistic  bottom  topography,  and 
coastline  geometry  that  is  determined  by  the  200-m  isobath. 
It  is  nonlinear,  thermodynamic,  eddy-resolving  (1/16°  res¬ 
olution  in  latitude  by  —0.088°  in  longitude)  and  has  a 
domain  extending  from  20°S  to  62°N  and  from 
109.125°E  to  77.203° W,  but  here  the  focus  is  on  the  GT 
region.  The  latitudinal  extent  of  the  model  domain  allows 
study  of  the  effects  of  equatorially  generated  signals  (i.e., 
CTWs)  in  the  eastern  Pacific  Ocean,  here  specifically  on  the 
GT  region.  The  model  is  forced  1979-2001  with  daily 
averaged  heat  fluxes  and  6-hourly  10-m  winds  from 
ECMWF  using  the  methodology  of  Kara  el  al.  [2002] 
and  does  not  include  ocean  data  assimilation. 

3.  Results  and  Discussion 

3.1.  Tehuantepec  Eddies  in  T/P  and  NLOM 

[fi]  The  T/P  satellite  altimeter  repeats  an  orbit  approxi¬ 
mately  every  10  days.  It  measures  sea  surface  height  (SSH) 
every  6.2  km  (1  second  average)  along  satellite  tracks.  The 
distance  between  two  parallel  tracks  is  about  300  km  in  our 
area  of  study.  This  wide  track  separation  usually  discour¬ 
ages  the  use  of  along-track  altimeter  measurements  to  track 
eddies.  However,  because  of  the  large  spatial  extent  of  the 
long-lived  Tehuantepec  eddies,  which  have  horizontal 
dimensions  of  more  than  300  km  (Figure  1),  and  westward 
propagation  speeds  of  —11-17  cm/s,  the  10-day  repeat 
sampling  of  T/P  depicts  their  passage  and  amplitude  very 
clearly  in  the  form  of  an  along-track  SSH  versus  time 
diagram  (positive  SSH  anomalies  in  Figure  2). 
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Figure  2.  (a)  Japanese  Meteorological  Agency  (JMA)  El  Nino  index.  Details  of  the  index  are  given  by 

JMA  [1991]  and  ftp://www.coaps.fsu.edu/pub/JMA  SST_lndex.  The  period  of  study  includes  the  three 
different  scenarios:  El  Nino  (1997),  La  Nina  (1998  and  1999),  and  neutral  years  (1992,  1993,  1994,  1995, 
1996,  and  2000).  Sea  surface  height  anomaly  time  series  from  (b)  TOPEX/Poseidon  (T/P)  and  (c)  1/16° 
Pacific  NLOM.  The  anomalies  are  relative  to  the  1993-1999  mean.  Units  are  centimeters  and  have  been 
calculated  along  the  bold  black  segment  of  the  T/P  satellite  track  shown  in  Figure  1.  (d)  Latitudinal 
average  of  the  time  series  in  Figures  2b  and  2c.  Blue  and  red  lines  correspond  to  T/P  and  NLOM, 
respectively.  A  1 -month  running  mean  filter  has  been  applied  to  the  latitudinally  averaged  time  series. 
The  correlation  coefficient  between  the  time  series  in  Figure  2d  is  0.72.  The  first  day  of  each  year  is 
labeled. 
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Figure  3.  Time  series  of  the  European  Centre  for  Medium-Range  Weather  Forecasts  (ECMWF)  10  m 
wind  stress  curl  (1  x  10~8  Pa/m)  root  mean  square  (RMS)  difference  with  respect  to  the  1979-2001 
mean  for  the  Gulf  of  Tehuantepec  and  Gulf  of  Papagayo  subregions  shown  in  Figures  6a  and  6d, 
respectively.  A  I -day  running  mean  has  been  applied  to  the  time  series.  The  RMS  for  each  of  the  eight 
time  series  is  indicated.  The  blue  and  red  time  series  correspond  to  the  Gulf  of  Tehuantepec  and 
Papagayo,  respectively.  The  First  day  of  each  year  is  labeled. 


[7]  The  T/P  and  NLOM  SSH  anomaly  data  for  the  period 
October  1992  to  March  2000  (Figures  2b  and  2c)  show  that 
three  to  five  anticyclonic  eddies  were  formed  during  each  of 
the  first  three  cold  seasons  (October- April)  (1992-1993, 
1993-1994,  and  1994-1995).  Particularly,  the  1992-1993 
anticyclonic  eddies  simulated  by  NLOM  were  weaker  than 
those  measured  by  T/P.  The  1995-1996  cold  season  was 
characterized  by  the  formation  of  three  anticyclonic  eddies. 
Throughout  the  fall  of  1996  no  eddies  were  generated; 
however,  in  the  1996-1997  winter  season  three  eddies 
formed.  The  T/P  and  NLOM  data  suggest  the  genesis  of 
at  least  seven  anticyclonic  Tehuantepec  eddies  during  the 
1997-1998  cold  season,  which  was  characterized  by  a 
strong  El  Nino  event  (Figure  2a).  Those  eddies  were 
characterized  by  the  largest  diameter  and  the  highest  SSH 
anomaly.  A  peculiar  characteristic  of  1997  is  that  the  eddy 
generation  continued  for  the  entire  year  in  both  T/P  and 
NLOM.  The  1998-1999  cold  season  T/P  data  record 
includes  just  two  anticyclonic  eddies  and  three  during  the 
1999-2000  cold  season.  During  these  last  two  cold  seasons 
a  La  Nina  event  was  taking  place  in  the  Pacific  Ocean 


(Figure  2a).  Table  1  summaries  the  number  of  Tehuantepec 
eddies  measured  by  T/P  and  simulated  by  NLOM  during  the 
8  cool  seasons  (October- April)  and  7  warm  seasons  (May- 
September)  analyzed  in  this  study. 

[s]  it  is  remarkable  how  the  model  (which  does  not 
include  ocean  data  assimilation)  reproduces  most  of  the 
Tehuantepec  eddy  characteristics  measured  by  T/P.  Indeed, 
the  correlation  coefficient  between  the  latitudinal ly  aver¬ 
aged  observed  (T/P)  and  simulated  (NLOM)  time  series  is 
0.72,  suggesting  a  semideterministic  nature  for  the  genera¬ 
tion  and  evolution  of  Tehuantepec  eddies  (Figure  2d). 

[9]  The  SSH  time  series  in  Figure  2  show  Tehuantepec 
eddy  formation  during  the  entire  year  of  1997,  even  when 
strong  Tehuantepec  wind  events  do  not  occur  (Figure  3). 
This  suggests  the  existence  of  an  alternative  generation 
mechanism  for  Tehuantepec  eddies,  and  advantages  of 
using  SSH  satellite  altimeter  data  instead  of  satellite  sea 
surface  temperature  (SST)  data  for  the  study  of  the  Tehuan¬ 
tepec  eddies.  CTWs  are  presented  as  the  alternate  mecha¬ 
nism,  a  topic  discussed  in  section  3.3.  The  warm-core 
anticyclonic  Tehuantepec  eddies  exist  in  a  region  where 
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Table  1.  Number  of  Anticyclonic  Tehuantepec  Eddies  Measured  by  T/P  and  Simulated  by  NLOM  During  the  Eight  Cool  Seasons 


(October-  April)  and  Seven  Warm  Seasons  (May-September)  Analyzed  in  This  Study3 


Cool  Season 

Warm  Season 

Total  in  T/P 

Total  in  NLOM 

Year 

Number 
of  Eddies 
in  T/P 

Number 
of  Eddies 
in  NLOM 

Year 

Number 
of  Eddies 
in  T/P 

Number 
of  Eddies 
in  NLOM 

1992  1993 

3 

3 

1993 

1 

1 

4 

4 

1993  1994 

4 

3 

1994 

l 

I 

5 

4 

1994  1995 

5 

4 

1995 

1 

0 

6 

4 

1995  1996 

3 

3 

1996 

0 

1 

3 

4 

1996  1997 

3 

3 

1997 

2 

2 

5 

5 

1997  1998 

7 

7 

1998 

1 

1 

8 

8 

1998  -1999 

2 

3 

1999 

2 

1 

4 

4 

1999  -2000 

3 

3 

:,The  total  number  of  anticyclonic  eddies,  in  T/P  and  NLOM,  includes  the  eddies  in  a  cool  season  plus  the  eddies  in  the  following  warm  season. 


the  surrounding  water  is  also  warm  [ Hansen  and  Maul , 
1991].  However,  during  their  formation  and  for  a  few  days 
after,  the  anticyclonic  Tehuantepec  eddies  are  clearly  rec¬ 
ognized  in  satellite  SST  images  (Figure  4).  This  is  enhanced 
by  the  strong  SST  gradient,  which  is  created  by  the  large 
mixing  and  upwelling  induced  SST  decrease  along  the  axis 
of  the  GT  wind  path.  Conditions  like  those  shown  in 
Figure  4  last  only  a  few  days  since  they  are  generated  by 
the  temporary  strong  wind  events.  When  the  eddies 
propagate  to  the  west  the  eddy  thermal  signature  is  difficult 
to  discern  in  satellite  SST  images  because,  once  the  strong 
cooling  ceases,  the  eddy  air-sea  temperature  differences 
stimulate  large  heat  fluxes  in  the  tropics  that  reduce  the 
strong  SST  gradient  and  eventually  render  the  eddies 
undetectable  in  thermal  images.  Therefore  satellite  altimeter 
SSH  observations  can  be  more  useful  for  the  study  of  eddies 
in  this  tropical  ocean. 

3.2.  GT  Wind  as  the  Generator  and  Modulator  of 
Tehuantepec  Eddies 

[to]  An  index  of  the  El  Nino  -  La  Nina  signal  (Figure  2a) 
and  the  SSH  anomaly  time  series  (Figures  2b-2d)  indicate 
that  the  seasons  characterized  by  the  fewest  and  the  smallest 
Tehuantepec  eddies  (1995-1996,  1998-1999,  and  1 999— 
2000)  coincided  with  the  cold  phase  of  the  index.  Whereas 
the  seasons  characterized  by  the  largest  (number,  diameter 
and  SSH  anomaly)  Tehuantepec  eddies  (1992-1993, 
1993-1994,  1994- 1 995,  and  1997- 1998)  occurred  during 
the  warm  phase.  One  explanation  for  this  modulation  of 
Tehuantepec  eddies  is  as  follows.  Given  the  results  in 
Figure  2  and  assuming  that  the  GT  winds  are  the  only 
mechanism  for  generation  and  modulation  of  Tehuantepec 
eddies,  then  fewer  and/or  weak  Tehuantepec  eddies  would 
be  caused  by  fewer  and/or  weaker  GT  wind  events.  This 
could  occur  as  a  result  of  weak  atmospheric  pressure 
gradients  between  the  Gulf  of  Mexico  and  the  Pacific  Ocean 
and  weaker  and/or  less  frequent  cold  fronts  passing  through 
the  Gulf  of  Mexico.  The  anomalous  location  of  the  jet 
stream,  which  is  shifted  southward  (northward)  during  the 
warm  (cold)  phase  of  the  index  strengthening  (weakening) 
the  cold  fronts  that  reach  the  Gulf  of  Mexico  [Quadrelli  and 
Wallace ,  2002],  is  consistent  with  Tehuantepec  eddy  inter¬ 
annual  variability  measured  by  T/P  (Figure  2b),  simulated 
by  NLOM  (Figure  2c)  and  with  the  preceding  explanation. 

[n]  To  test  this  hypothesis  we  analyzed  1979-2001 
ECMWF  6-hourly  10-m  winds,  1979-2001  ECMWF  12- 
hourly  1000  mb  winds  and  1990-2001  NOGAPS  6-hourly 


surface  stresses  over  the  GT  region.  An  example  for  the 
ECMWF  10-m  winds  over  the  period  October  1992  to 
September  2000  is  shown  in  Figure  3.  The  analyses  show 
that  the  ECMWF  and  NOGAPS  winds  reproduce  several 
satellite  scatterometer  and  in  situ  measured  features  of  the 
GT  winds  [B  our  ass  a  el  a/.,  1999;  Chelton  et  al ,  2000a, 
2000b;  Romero-Centeno  et  al .,  2003]  (http://ww\v.coaps. 
fsu.edu/qscat/anim).  Those  characteristics  are  strong  wind 
events  mainly  during  fall-winter,  a  wind  event  duration  of 
—2-3  days,  and  wind  speed  and  negative  (positive)  wind 
stress  curl  on  the  west  (east)  side  of  the  wind  path. 
Furthermore,  Kelly  et  al.  [1999]  performed  an  exhaustive 
model/data  comparison  between  ECMWF  winds  and  the 
NASA  scatterometer  (NSCAT)  winds  in  the  tropics,  includ¬ 
ing  the  Tehuantepec  and  Central  America  regions,  which 
shows  the  capability  of  ECMWF  to  accurately  represent  the 
winds  of  the  region.  Even  though  the  different  wind 
products  realistically  represent  the  GT  winds,  the  statistics 
(the  root  mean  square  of  the  wind  stress  curl)  do  not  show 


: 


Figure  4.  Sea  surface  temperature  (°C)  in  the  Gulf  of 
Tehuantepec  region  as  measured  by  the  NOAA  14  satellite 
on  22  January  1996.  Image  processed  by  Agustin  Fernandez 
(Universidad  Nacional  Autonoma  de  Mexico).  It  shows 
upwelling  and  mixing  (cold  SST)  and  an  anticyclonic  eddy 
driven  by  a  Tehuantepec  wind  event. 
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Figure  5.  Sea  surface  height  anomaly  (color  contours  in  centimeters)  for  nine  different  dates  in 
February- March  1997  as  simulated  by  1/16°  Pacific  NLOM  forced  by  ECMWF  6  hourly  10  m  winds 
and  daily  averages  for  the  other  fields  used  in  the  thermal  forcing.  The  locations  of  the  Gulf  of  Papagayo 
(GP),  Gulf  of  Tehuantepec  (GT),  and  Acapulco  (ACA)  are  indicated.  The  gray  area  represents  ocean 
areas  not  modeled  by  Pacific  NLOM,  i.e.,  depths  shallower  than  200  m. 


significant  interannual  variability  (Figure  3).  In  addition,  a 
visual  inspection  of  the  time  series  in  Figure  3  does  not 
indicate  important  differences  in  the  number  and  intensity  of 
the  wind  events  during  the  two  opposite  cool  season 
regimes,  the  1997-1998  (El  Nino  year)  and  1999-2000 
(La  Nina  year).  Note  also  that  during  1997  the  Tehuantepec 
eddy  generation  extends  beyond  the  1996-1997  cold  sea¬ 
son.  That  fact  cannot  be  explained  in  terms  of  wind  stress 
curl  alone,  because  the  wind  stress  curl  in  spring-summer  is 
significantly  weaker  than  in  fall-winter  (Figure  3).  Hence 
the  lack  of  interannual  variability  in  the  wind  fields  suggests 
that  by  themselves  the  GT  winds  are  not  sufficient  to 
explain  the  satellite  measured  Tehuantepec  eddy  interarmual 
variability. 

3.3.  Coastally  Trapped  Waves  and  the  Tehuantepec 
Eddies 

3.3.1.  Coastally  Trapped  Waves  and  Winds 

[12]  A  sequence  of  snapshots  showing  the  GT  winds  and 
the  CTW  as  modulators  and  generators  of  Tehuantepec 
eddies  are  presented  in  Figures  5  and  6.  On  15  February 
1997  the  SSH  anomaly  around  the  GT  is  characterized  by 
two  anticyclonic  eddies  and  one  cyclonic  eddy  (Figure  5a). 


The  anticyclones  have  maximum  SSH  anomalies  of  ~15 
and  20  cm  and  centers  close  to  95°W,  1  L8°N  and  97.5°W, 
14°N,  respectively.  The  cyclonic  eddy  has  minimum  SSH 
anomaly  of  11  cm  at  its  center  (94. 5° W,  13.6°N).  The 
three  eddies  were  generated  by  the  wind  events  that  oc¬ 
curred  in  the  GT  during  January  1997  (Figure  3).  In 
addition,  Figure  5a  shows  areas  of  positive  SSH  anomaly 
around  90°W,  12°N.  The  wind  stress  curl  is  strong  and 
anticyclonic  (cyclonic)  on  the  western  (eastern)  side  of  the 
GT  (Figure  6a).  This  wind  stress  curl  strengthened  the  larger 
anticyclonic  eddy  and  the  cyclonic  eddy,  increasing  the 
amplitude  of  the  SSH  anomalies  to  a  maximum  >25  cm  and 
a  minimum  <—15  cm,  respectively  (Figures  5b).  In  addi¬ 
tion,  the  tip  of  a  downwelling  CTW  is  entering  the  region 
near  85°W,  10°N  and  there  is  weak  wind  stress  curl  in  the 
entire  study  area,  except  for  the  GT,  as  can  be  seen  in 
Figures  5b  and  6b.  Next,  the  wind  stress  curl  weakened 
considerably,  the  CTW  propagated  northward,  and  a  new 
anticyclonic  eddy  formed  on  the  southeast  side  of  the  GT 
(92.5°W,  I2.5°N)  (Figures  5c-5g  and  6c-6g).  Note  that  in 
spite  of  the  weak  anticyclonic  wind  stress  curl  to  the 
southeast  of  the  GT,  which  is  ~7  times  smaller  than  to 
the  west  of  the  GT  (Figure  6),  an  anticyclonic  eddy  was 
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Figure  6.  Wind  stress  curl  (color  contours  in  1  x  1(F8  Pa/m)  for  nine  different  dates  in  February - 
March  1997  as  determined  from  the  ECMWF  10-m  winds.  The  magenta  boxes  in  Figures  6a  and  6d 
indicate  the  two  areas  where  the  time  series  of  Figure  3  were  extracted. 


generated  in  this  region  (Figure  5).  Finally,  the  tail  of  the 
CTW  reached  and  crossed  the  Acapulco  region  and  the  area 
of  study  was  populated  with  larger  and  stronger  eddies  than 
a  month  prior  (Figures  5h-5i). 

[13]  It  is  noteworthy  that  the  GT  winds  strengthened  the 
anticyclonic  (cyclonic)  eddy  located  on  the  western  (central- 
eastern)  side  of  the  GT  on  1 5  February  1 997  (Figures  5  and  6). 
That  is  an  example  of  the  GT  winds  as  a  generator  and 
modulator  of  the  Tehuantepec  eddies.  However,  equally 
notable  is  the  formation  of  the  anticyclonic  eddy  to  the 
southeast  of  the  GT  that  coincides  with  the  passage  of  the 
CTW  through  the  area  during  a  period  of  weak  anticyclonic 
wind  stress  curl  (Figures  5c-5g  and  6c-6g).  Hence  some 
questions  present  themselves.  What  is  the  origin  of  these 
downwelling  CTWs?  How  do  they  fluctuate  interannually? 
How  do  the  CTWs  contribute  to  the  modulation/formation  of 
the  Tehuantepec  eddies? 

[m]  During  the  boreal  fall/winter  season  the  equatorial 
Pacific  is  characterized  by  the  formation  of  2-4  intra- 
seasonal  (30-90  day)  baroclinic  downwelling  Kelvin 
waves,  which  are  more  (less)  common  and  more  intense 
(weaker)  during  El  Nino  (La  Nina)  years.  These  waves 
propagate  eastward  until  they  reach  the  west  coast  of  the 
Americas,  where  they  split  into  northward  and  southward 


propagating  CTWs  with  the  northward  CTWs  reaching  the 
GT  and  beyond  [Spiilane  el  al.,  1987;  Enfield ,  1987; 
Kessler  et  al. ,  1995;  Melsom  et  al. ,  2003].  Figure  7  is  a 
diagram  of  along  equatorial  and  along  coastal  SSH  versus 
time,  as  simulated  by  NLOM.  It  shows  (1)  the  western  and 
central  Pacific  as  the  generation  region  for  the  equatorial 
waves,  (2)  a  nearly  constant  propagation  speed  (~2.3  m/s) 
for  the  waves,  which  denotes  their  nondispersive  character, 
(3)  interannual  variability  in  both  number  and  strength  of 
the  waves,  (4)  proliferation  of  the  waves  during  the  boreal 
cold  season  and  during  El  Nino  events,  and  (5)  wave 
amplification  as  the  waves  propagate  poleward  along  the 
coast. 

[15]  Considering  that  both  the  intraseasonal  baroclinic 
downwelling  CTWs  and  the  Tehuantepec  winds  are  distinc¬ 
tive  features  of  the  eastern  tropical  Pacific,  which  are 
maximal  during  the  boreal  cool  season,  the  oceanic  effects 
of  both  may  coincide  in  space  and  time  around  the  GT.  The 
Tehuantepec  winds,  which  inject  large  amounts  of  energy 
into  the  GT  waters  (Figure  8a),  induce  negative  (positive) 
relative  vorticity  on  the  western  (central-eastern)  side  of  the 
GT  (Figure  8b)  via  Ekman  pumping  (suction).  In  the 
process  a  strong  converging  (diverging)  Ekman  transport 
increases  (decreases)  the  oceanic  pressure  (Figure  8c)  and  in 
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Figure  7.  Sea  surface  height  anomaly  (color  contours  in  centimeters)  time  series  from  1/16°  Pacific 
NLOM,  first  along  the  equator  (starting  in  the  western  Pacific  and  propagating  eastward  until  arrival  at 
the  Americas  West  Coast),  and  second  along  the  coast  to  the  20°N.  The  white  vertical  lines  to  the  west  of 
130°E  indicate  the  location  of  Equatorial  Western  Pacific  Islands  and  the  white  vertical  line  close  to 
270°E  indicates  the  location  of  the  Galapagos  Islands.  The  equatorial  position  of  the  Americas  West 
Coast  is  indicated  by  AME.  The  first  day  of  each  year  is  labeled. 


the  convergence  (divergence)  areas,  lowers  (raises)  the 
thermocline  (Figure  8d),  increases  (decreases)  the  SSH 
(Figure  8e),  decreases  (increases)  potential  vorticity 
(Figure  80,  and  forms  anticyclonic  (cyclonic)  eddies 


(Figure  8).  At  the  same  time,  the  northward  propagating 
baroclinic  downwelling  CTW  induces  strong  poleward 
currents  as  it  passes  (Figure  9),  increasing  the  energy  along 
the  coast  (Figure  8a).  Also,  this  wave  increases  the  coastal 
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Figure  8.  Snapshots  for  1  January  1997  as  simulated  by  1/16°  Pacific  NLOM.  (a)  Surface  kinetic 
energy  anomaly  per  mass  (color  contours  in  log  (m2/s2)),  (b)  surface  relative  vorticity  anomaly  (color 
contours  in  I  x  I0~5  s  !),  (c)  density  normalized  surface  pressure  anomaly  (color  contours  in  m2/s2) 
and  currents  (arrow  vectors),  (d)  surface  layer  thickness  (color  contours  in  meters),  (e)  sea  surface 
height  (color  contour  in  centimeters),  and  (f)  surface  potential  vorticity  anomaly  (color  contours  in  1  x 
10  6  m  1  s  l). 
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Figure  9.  Surface  currents  (arrow  vectors)  and  current  speed  (color  contours  in  cm/s)  for  1  January 
1997  as  simulated  by  1/16°  Pacific  NLOM. 


SSH  (Figure  8e),  increases  the  coastal  pressure  (Figure  8c), 
and  depresses  the  coastal  thermocline  (Figure  8d).  Further¬ 
more,  the  currents  induced  by  the  CTW  produce  a  strong 
horizontal  shear  flow  that  is  a  source  of  negative  (positive) 
relative  vorticity  in  the  inshore  (offshore)  oceanic  regions 
during  the  passage  of  the  wave  (Figures  8b  and  9).  This 
extra  negative  relative  vorticity  is  restricted  to  the  lateral 
frictional  boundary  layer,  which  has  an  offshore  extension 
of  MO  km  and  is  several  times  smaller  than  the  radius  of 
the  anticyclonic  eddies  (>100  km)  (Figures  8  and  9).  Hence 
this  extra  negative  relative  vorticity  is  probably  not  respon¬ 
sible  for  initiating  the  formation  of  the  eddies.  In  addition, 
the  northward  propagating  downwelling  CTW  can  be  a 
source  of  low  values  of  potential  vorticity  for  the  GT  region. 
As  an  example,  note  the  similarity  between  the  potential 
vorticity  values  in  the  equatorial  Pacific  and  the  potential 
vorticity  values  along  the  coast,  which  are  associated 
with  the  passage  of  the  CTW  through  the  region,  and  the 
potential  vorticity  values  of  the  anticyclonic  eddies 
(Figure  10).  How  do  intraseasonal  baroclinic  CTWs  influ¬ 
ence  on  the  generation  of  the  eddies?  This  question  is 
addressed  in  the  next  section. 

3.3.2.  Intraseasonal  CTWs  and  Eddy  Generation 

[i6]  According  to  linear  equatorial  wave  theory,  for  a 
given  frequency  there  exists  a  critical  latitude  below  which 
unforced  CTWs  radiate  energy  offshore  in  the  form  of  long 
Rossby  waves  [Clarke  and  Shi,  1991].  Considering  that  for 


intraseasonal  frequencies  (30-90  days)  the  nonmeridional 
coasts  of  Central  America  and  southwest  Mexico  are 
located  to  the  north  of  the  critical  latitude  [Clarke  and 
Shi ,  1991,  Table  1],  then  the  free  propagating  intraseasonal 
CTWs  of  Figures  5,  7,  and  8  will  not  separate  from  the 
coast.  Consequently  the  formation  of  eddies  via  CTW 
detachment  from  the  coasts  of  Central  America  and  south¬ 
west  Mexico  is  not  a  mechanism  that  will  occur  during  the 
passage  of  intraseasonal  CTWs.  Also,  for  the  particular  case 
of  the  anticyclonic  eddy  centered  near  92.5°W,  13.2°N  in 
Figure  5f,  the  wind  as  the  eddy  generator  can  be  dismissed, 
since  this  eddy  formed  during  a  period  of  weak  winds 
(Figure  6).  Then,  what  was  the  generation  mechanism  for 
this  eddy? 

[i?]  The  CTWs  discussed  here  are  generated  by  equato¬ 
rial  winds.  Thus  they  could  be  surface  intensified 
and  characterized  by  a  vertical  shear  flow  (Figures  5,  9, 
and  11),  making  them  candidates  for  triggering  baroclinic 
instability.  Also,  CTWs  induce  strong  surface  currents  in  the 
areas  of  their  passage,  generating  a  strong  horizontal  shear 
flow  (Figures  5,  9,  and  11),  making  them  a  candidate  for 
triggering  barotropic  instability.  Thus  we  hypothesize  that 
Tehuantepec  eddies  can  be  modulated  and/or  generated  by 
barotropic  and/or  baroclinic  instabilities,  which  are  trig¬ 
gered  by  CTWs. 

[is]  Baroclinic  instabilities  occur  when  perturbations 
grow  by  drawing  energy  from  the  available  potential  energy 
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Figure  10.  Surface  potential  vorticity  anomaly  (color  contours  in  1  x  10  6  m  1  s  ')  for  1  January  1997 
as  simulated  by  1/16°  Pacific  NLOM. 


associated  with  sloping  isopycnals  [, Kundu ,  1 990].  In  the  case 
of  baroclinic  downwelling  CTWs  these  are  characterized  by  a 
sharp  rise  of  the  SSH  toward  the  coast  (Figures  5b-5i). 
Spall  [1992]  proposed  and  calculated  baroclinic  instabilities 
in  the  southwestward  eastern  boundary  current  known  as 
the  Cape  Verde  Frontal  Zone  as  a  means  to  generate 
radiating  baroclinic  Rossby  waves,  which  produce  meso- 
scale  variability  in  the  mid-ocean.  In  the  GT  case  the 
current  is  a  northwestward  eastern  boundary  current,  which 
is  associated  with  the  passage  of  poleward  propagating 
baroclinic  downwelling  CTWs  that  could  trigger  the  devel¬ 
opment  of  barotropic  and  baroclinic  instabilities,  which  in 
turn  could  trigger  the  development  of  eddies  and  related 
mesoscale  variability.  Thus  we  search  for  signatures  of 
barotropic  and  baroclinic  instabilities  in  the  pressure,  energy 
(section  3.3.3)  and  currents  fields,  and  we  isolate  the  effects 
of  baroclinic  instabilities  analyzing  the  results  of  a  1.5-layer 
reduced  gravity  simulation  (section  3.3.3). 

[d]  Some  support  for  the  hypothesis  of  baroclinic  insta¬ 
bilities  as  a  trigger  for  development  of  Tehuantepec  eddies 
is  provided  by  the  vertical  structure  of  the  eddies.  As  an 
example  of  that,  the  surface  layer  pressure  overlain  on  the 
fifth  layer  pressure  (Figure  12)  shows  a  significant  phase 
difference  between  the  two  layers  in  the  vicinity  of  the  eddy 
centered  near  93°W,  13.2°N.  This  is  a  signature  of  baro¬ 
clinic  instability  when  it  occurs  during  a  period  of  rapid 
eddy  intensification,  but  can  occur  for  any  propagating 


baroclinic  eddy  because  of  the  pattern  of  vortex  stretching 
and  compression  [e.g.,  see  Hurlburt  et  al. ,  1990]. 

[20]  What  is  the  role  of  barotropic  instabilities  in  the 
development  of  the  Tehuantepec  eddies?  Was  the  new  eddy 
in  Figure  5  triggered  by  barotropic  instabilities?  Barotropic 
instabilities  are  due  to  perturbations  where  the  energy  is 
supplied  by  the  horizontal  shear  of  the  currents.  To  provide 
some  insight  into  the  origin  of  the  instabilities  generated  by 
the  CTW  in  Figure  5,  a  Beta  Rossby  number,  which  is 
defined  as  the  ratio  of  relative  to  planetary  vorticity  advec- 
tion,  Rft  =  'o/fSr2,  was  calculated.  Here  v  is  the  maximum 
swirl  velocity  of  an  eddy  averaged  around  the  eddy,  and  r  is 
the  mean  radius  of  the  eddy  alv.  Rs  of  order  1(10)  suggests 
barotropic  (baroclinic)  instabilities  involving  the  first  bar¬ 
oclinic  (barotropic)  mode  [McWilliams  and  Flierl ,  1979; 
Hurlburt  and  Thompson ,  1982,  1984;  Murphy  et  al .,  1999]. 
In  the  case  of  the  eddy  with  center  close  to  93°W,  1 3.2°N  on 
12  March  1997  (Figure  5h)  0  ^  2.2  x  10"  11  rrT1  s'1,  *u  = 
50  cm/s  and  r  ^  100  km.  Consequently,  RB  ^  2.3  suggest¬ 
ing  a  key  role  for  barotropic  or  internal  mode  baroclinic 
instabilities  [Hurlburt  et  al. ,  1996]  in  the  formation  of  this 
particular  eddy. 

[21]  Taken  together,  the  results  in  Figures  5,  6,  and  11-12 
show  signatures  consistent  with  barotropic  and  baroclinic 
instabilities  triggered  by  a  short-wavelength  intraseasonal 
baroclinic  downwelling  CTW  during  a  period  of  weak  GT 
winds,  and  the  possible  influence  of  the  instabilities  on  the 
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Figure  11.  Currents  (arrow  vectors)  and  current  speed  (color  contours  in  cm/s)  for  six  layers  on  3  March 
1997  as  determined  by  1/16°  Pacific  NLOM. 
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Figure  12.  Surface  layer  pressure  (red  in  m2/s2)  overlain  on  layer  5  pressure  (black  in  m2/s2)  for 
12  March  1997  as  simulated  by  1/16°  Pacific  NLOM.  Note  the  phase  lag  in  the  vicinity  of  the  eddy 
centered  near  93°W,  13.2°N. 
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development  of  an  anticyclonic  eddy  southeast  of  the  GT. 
Since  this  new  mechanism  for  generation  of  Tehuantepec 
eddies  is  not  wind  dependent,  then  it  should  be  possible  to 
isolate  its  effects  during  summer  time,  when  strong  offshore 
GT  wind  events  are  rare  (Figure  3). 

3.3.3.  Tehuantepec  Eddy  Generation  During  the 
Summer  of  1997 

[22]  To  provide  additional  insight  concerning  CTW  im¬ 
pact  on  Tehuantepec  eddy  formation,  we  document  the 
generation  of  two  anticyclonic  eddies  during  the  summer 
of  1997.  These  occur  during  a  5-month  period  (mid-May  to 
early  October)  of  weak  GT  and  Papagayo  winds  (Figure  3). 
Detailed  discussion  of  the  Papagayo  winds  is  given  by 
Schultz  et  al.  [1997,  and  references  therein]  and  Chelton 
et  al.  [2000a,  2000b,  and  references  therein].  On  14  June 
1997  the  SSH  anomaly  field  along  the  coast  from  Central 
America  to  the  entrance  of  the  Gulf  of  California  was 
characterized  by  generally  weak  highs  and  lows  in  SSH, 
generally  weak  cyclonic  and  anticyclonic  eddies,  a  lack  of 
anticyclonic  eddies  in  the  GT  area,  and  by  only  a  weak 
CTW  (Figure  13a).  However,  two  weeks  later  the  tip  of  an 
arriving  downwelling  CTW  (positive  SSH  anomaly  of 
~20  cm)  can  be  seen  along  the  Central  America  coast 
(85°W,  I0°N)  in  Figure  1 3b.  This  long-wavelength  (i.e., 
south  of  the  critical  latitude)  CTW  (which  is  associated  with 
the  1997  El  Nino  event)  continued  propagating  northward 
and  raising  the  coastal  region  SSH  from  Central  America  to 
the  entrance  of  the  Gulf  of  California  (Figures  13b-13g). 
According  to  the  linear  equatorial  wave  theory  for  interan¬ 


nual  frequencies,  the  nonmeridional  coasts  of  Central 
America  and  southwest  Mexico  are  located  to  the  south  of 
the  critical  latitude  [Clarke  and  Shi,  1991,  Table  I],  Hence, 
during  its  northward  propagation  the  interannual  CTW 
radiates  energy  westward  in  the  form  of  Rossby  waves 
(Figure  13).  The  simultaneous  occurrence  of  these  two 
processes,  the  northward  propagation  of  the  CTW  and  its 
westward  radiation  of  energy,  can  be  recognized  by  noting 
the  narrow  cross-shore  trapping  scale  of  the  CTW  (which  is 
proportional  to  the  radius  of  deformation)  (Figures  13b  — 
13g)  and  the  subsequent  westward  propagation  of  this 
positive  SSH  anomaly.  By  17  August  1997  it  is  completely 
separated  from  the  coast  (Figures  I3h-13r).  Also,  note  that 
on  1 8  July  1997  a  short-wavelength  baroclinic  downwelling 
CTW  (which  is  characterized  by  a  SSH  anomaly  of  ~30  cm) 
reached  the  area  of  study  and  its  passage  through  the  GT 
region  coincided  with  the  formation  of  anticyclonic  Tehuan¬ 
tepec  eddies  on  the  southeast  (90° W,  12.7°N)  and  northwest 
(96.5°W,  15.5°N)  sides  of  the  GT  (Figure  13h-13n).  Since 
these  two  eddies  were  formed  during  the  5-month  period  of 
weak  GT  and  Papagayo  winds  during  1997  (Figure  3),  the 
wind  can  be  dismissed  as  the  source  of  energy  for  the 
generation  of  the  eddies.  Then,  what  is  the  generation 
mechanism  for  these  eddies?  Did  the  CTWs  in  Figure  13 
trigger  barotropic  and/or  baroclinic  instabilities? 

[23]  The  propagation  of  the  two  CTWs  along  the  coasts 
of  Central  America  and  southwest  Mexico  is  illustrated  by 
the  thickness  of  the  two  top  model  layers  (Figure  14).  The 
coastal  deepening  effects  of  the  long-wavelength  CTW  are 
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Figure  14.  Layer  1  plus  layer  2  thickness  (color  contours  in  m)  for  six  different  dates  in  July  1997  as 
simulated  by  1/16°  Pacific  NLOM. 


Figure  13.  SSH  anomaly  (color  contours  in  cm)  for  18  different  dates  in  June-August  1997  as  simulated  by  1/16°  Pacific 
NLOM.  The  positions  of  the  Gulf  of  Tehuantepec  (GT),  Acapulco  (ACA),  Cabo  Corrientes  (CC),  and  Maria  Islands  (Ml) 
are  indicated.  The  magenta  box  in  Figure  1 3 h  indicates  the  area  where  the  kinetic  energy  time  series  of  Figure  15  were 
calculated. 
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Figure  15.  Time  series  of  kinetic  energy  averaged  over  the 
magenta  box  shown  in  Figure  13h.  The  blue,  red,  and  green 
time  series  correspond  to  the  second  layer  (IQ)  and  the  two 
lowest  layers  (IQ,  and  IQ),  respectively. 

clearly  depicted  in  Figures  14a- I4f.  On  15  July  1997  the 
top  two  layer  thickness  includes  a  narrow  and  elongated 
feature  (which  is  associated  with  the  long-wavelength 
CTW)  that  is  maximum  along  the  coast,  decays  offshore, 
and  is  characterized  by  cross-shore  and  alongshore  exten¬ 
sions  of  ^80  and  ~3000  km  for  the  positive  anomaly, 
respectively  (Figure  14a).  Also,  there  is  an  anticyclonic 
eddy  offshore  centered  near  103°W,  12°N.  However,  there 
is  no  anticyclonic  eddy  activity  in  the  coastal  areas  north¬ 
west  and  southeast  of  the  GT.  Next,  the  short-wavelength 
CTW  propagates  through  the  GT  region  (Figures  14b-  14e). 
Coincident  with  the  passage  of  a  short-wavelength  CTW 


through  the  GT  region,  the  long-wavelength  CTW  starts  to 
meander  (Figure  14c).  It  continues  in  an  amplifying  mean¬ 
der  state  (Figures  14c- 14f)  until  anticyclonic  eddies  are 
generated  northwest  and  southeast  of  the  GT  (Figure  14f). 
The  alongshore  wavelength  of  the  meander  is  ~700  km 
(Figure  I4e).  That  wavelength  could  be  indicative  of 
baroclinic  instability.  Kundu  [1990]  shows  that  in  a  2-layer 
system,  the  wavelength  with  the  fastest  growth  rate  is  ~3.9  • 
i \  •  R\  (where  R]  is  the  first  baroclinic  radius  of  deforma¬ 
tion),  which  is  ~735  km  in  our  area  of  study,  similar  to  the 
~700  km  wavelength  simulated  by  NLOM  (Figure  14e). 

[24]  As  it  propagates  away  from  the  coast,  the  long- 
wavelength  CTW  develops  meanders  at  a  wavelength 
consistent  with  baroclinic  instability  involving  the  baro- 
tropic  mode.  However,  the  evidence  points  toward  internal 
mode  baroclinic  instability  for  the  actual  growth  of  the 
eddies,  occurring  when  the  short  CTW  passes  through  the 
area.  This  evidence  is  provided  by  kinetic  energy  and 
pressure  fields.  Figure  15  shows  time  series  of  kinetic 
energy  in  the  second  layer  (K2)  and  the  lowest  two  layers 
(IQ,  and  IQ)  during  July  1997  for  a  region  in  the  vicinity  of 
the  eddy  southeast  of  the  GT  as  outlined  in  Figure  13h.  The 
large  K2  peak  on  18  July  is  associated  with  the  passage  of 
the  short-wavelength  CTW  through  the  region.  It  is  note¬ 
worthy  that  of  all  the  peaks  in  K2,  only  the  peak  on  18  July 
has  a  corresponding  increase  in  IQ,  which  is  not  reflected  in 
IQ  (Figure  15).  The  decreasing  K2  accompanied  by  increas¬ 
ing  IQ  and  a  basically  invariant  IQ  indicates  a  rapid  transfer 
of  energy  from  the  upper  to  the  lower  baroclinic  layers,  a 
signature  typical  of  internal  mode  baroclinic  instability  as 
clearly  demonstrated,  for  example,  by  Holland  and  Lin 
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Figure  16.  Surface  layer  pressure  anomaly  (color  contours  in  m2/s2)  overlain  on  layer  5  pressure  (black 
(negative)  and  white  (positive)  contours  in  m2/s2)  for  30  July  1997  as  simulated  by  1/16°  Pacific  NLOM. 
Note  the  phase  lag  in  the  vicinity  of  the  eddies  centered  near  97°W,  14.8°N  and  91.8°W,  12.0°N. 
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[1975]  and  Hurlburt  et  al.  [1990]  for  baroclinic  instability 
involving  the  barotropic  mode.  The  increase  in  K5  on 
18  July  is  a  not  a  result  of  energy  entering  the  subdomain 
through  the  boundaries.  Furthermore,  the  surface  layer 
pressure  overlain  on  the  fifth  layer  pressure  (Figure  16) 
shows  a  significant  phase  difference  between  the  two  layers 
in  the  vicinity  of  the  eddies  centered  near  97°  W,  14.8°N, 
and  91.8°W,  12.5°N,  again  a  characteristic  of  baroclinic 
instability.  However,  because  of  the  pattern  of  vortex 
stretching  and  compression,  propagating  baroclinic  eddies 
in  general  show  a  phase  shift  with  depth.  Thus  we  also 
investigate  whether  barotropic  instability  could  be  active  in 
the  formation  of  these  eddies. 

[25]  The  Beta  Rossby  numbers  of  the  eddies  centered 
near  91.8°W,  12.5°N  (v  =  100  cm/s  and  r  ^  146  km),  and 
97°W,  I4.8°N  (v  =  60  cm/s  and  r  ^  160  km)  are  2.0,  and 
1 .2,  respectively.  Those  numbers  suggest  a  possible  role  for 
barotropic  instabilities  in  the  generation  of  these  two  eddies 
[Hurlburt  and  Thompson ,  1982],  but  Hurlburt  et  al.  [1996] 
found  similar  values  in  an  eddy-resolving  5.5-layer  reduced 
gravity  model  which  allows  internal  mode  baroclinic  insta¬ 
bility  but  excludes  the  barotropic  mode.  Then,  to  investigate 
further  the  possibility  of  barotropic  instability,  we  per¬ 
formed  a  simulation  using  a  1.5-layer  reduced  gravity 
configuration.  One  of  the  big  advantages  of  this  simulation 
is  that  it  does  not  include  baroclinic  instabilities  but  includes 
barotropic  instabilities,  allowing  us  to  isolate  the  effects  of 
the  barotropic  instabilities  in  the  Tehuantepec  eddy  gener¬ 
ation  and  modulation.  The  results  of  the  1.5-layer  reduced 
gravity  simulation  show  the  formation  of  the  two  anticy- 
clonic  Tehuantepec  eddies  during  July  1997  (Figure  17). 
Those  eddies  are  weaker  than  the  corresponding  eddies  in 
the  7-layer  simulation  with  the  maximum  swirl  velocity 
reduced  to  about  1/2.  The  Beta  Rossby  numbers  of  the 
eddies  centered  near  93° W,  12°N  (v  =  60  cm/s  and  r  zz 
1 50  km),  and  97.5°W,  1 5°N  (u  =  30  cm/s  and  r  ^  125  km) 
in  Figure  17  are  1.2,  and  1.0,  respectively.  Why  are  the 
eddies  weaker  in  the  1.5-layer  reduced  gravity  simulation 
when  the  eddies  cannot  pump  energy  to  any  other  vertical 
mode  as  in  the  7-layer  simulation,  and  when  consequently 
one  might  expect  the  generation  of  stronger  rather  than 
weaker  eddies  in  the  1.5-layer  reduced  gravity  simulation? 
The  eddies  in  the  1.5-layer  reduced  gravity  simulation  are 
weak  because  they  were  not  amplified  by  baroclinic  insta¬ 
bilities,  like  the  ones  which  imprint  their  signature  in  the 
kinetic  energy  time  series  of  Figure  15.  However,  because 
these  eddies  developed,  barotropic  instabilities  may 
contribute  to  the  development  of  the  anticyclonic  eddies 
centered  near91.8°W,  12.5°N  and  97° W,  14.8°N  on  30  July 
1997  (Figure  13)  as  well  as  baroclinic  instabilities. 
Therefore,  in  this  new  generation  mechanism  for  Tehuante¬ 
pec  eddies,  the  model  results  show  that  CTWs  can  trigger 
the  development  of  mixed  barotropic- internal  mode- 
baroclinic  instabilities  and  the  instabilities  trigger  the 
development  and/or  amplification  of  the  eddies. 


[26]  The  contributions  of  long-wavelength  and  short- 
wavelength  CTWs  to  the  formation  of  Tehuantepec  eddies 
during  the  summer  of  1997  are  summarized  as  follows. 
First,  the  long-wavelength  CTW  sets  up  the  background 
conditions  by  deepening  the  coastal  thermocline  (Figure  14) 
and  inducing  poleward  coastal  currents.  Next,  the  short- 
wavelength  CTW  strengthens  the  coastal  poleward  currents, 
producing  strong  horizontal  and  vertical  shear  of  the  veloc¬ 
ity  and  triggering  barotropic  and  baroclinic  instabilities 
(Figures  15  and  16).  This  causes  the  long-wavelength 
CTW  meanders  to  break,  generating  the  two  eddies 
(Figures  13  and  14).  Since  the  long-wavelength  CTW 
radiates  energy  westward  in  the  form  of  Rossby  waves,  it 
contributes  to  the  rapid  detachment  of  the  eddies  from  the 
coast,  aided  by  an  upwelling  CTW  which  arrives  around 
14  August  (Figure  13).  The  model  results  on  CTWs  as  gen¬ 
erators  of  Tehuantepec  eddies  are  corroborated  by  model- 
independent  analyses  of  satellite  altimeter  SSH  data.  The 
SSH  anomalies  from  the  Modular  Ocean  Data  Assimilation 
System  [Fox  et  al .,  2002]  clearly  show  the  arrival  of  the 
interannual  CTW  and  the  generation  of  two  eddies  southeast 
of  the  GT  during  summer  of  1 997  (Figure  1 8),  results  similar  to 
the  SSH  anomalies  simulated  by  NLOM  (Figure  13). 

3.3.4.  Interannual  Variability  of  Tehuantepec  Eddies 

[27]  The  GT  winds  as  the  main  generator  of  the  Tehuan¬ 
tepec  eddies  is  unquestionable.  On  the  other  hand,  the 
Tehuantepec  and  Papagayo  winds  do  not  show  significant 
interannual  variability  (Figure  3),  while  the  observed  and 
modeled  Tehuantepec  eddies  (Figure  2  and  Table  1)  and  the 
intraseasonal  CTWs  do  [ Enfield ,  1987;  Kessler  et  al 1995]. 
Downwelling  CTWs  contribute  to  the  generation  and/or 
amplification  of  eddies  by  creating  a  more  favorable  ambi¬ 
ent  environment  (deeper  thermocline  and  larger  positive 
SSH  anomalies)  and  by  triggering  baroclinic  and/or  baro¬ 
tropic  instabilities  (sections  3.3.2  and  3.3.3).  Hence  these 
new  generation  mechanisms  for  Tehuantepec  eddies  have  an 
intrinsic  interannual  variability  (Table  1),  which  has  been 
measured  indirectly  by  the  T/P  altimeter  and  directly 
simulated  by  NLOM  (Figure  2). 

[28]  El  Nino  and  La  Nina  events,  documented  in  Figure  2, 
are  two  opposite  examples  that  indicate  the  importance  of 
CTWs  as  a  modulator  of  Tehuantepec  eddy  development. 
While  the  number  and  intensity  of  anticyclonic  Tehuantepec 
eddies  and  intraseasonal  downwelling  CTWs  increases 
(decreases)  during  El  Nino  (La  Nina)  years  [Kessler  et  al 
1995]  (Figure  2  and  Table  1),  the  GT  and  Papagayo  winds 
are  not  significantly  different  during  those  events  (Figure  3). 
Also,  the  mesoscale  features  documented  in  Figures  5-18 
are  examples  that  illustrate  the  roles  of  CTWs  in  Tehuante¬ 
pec  eddy  formation  and/or  amplification.  However,  these 
are  not  unusual  examples.  Several  other  modeled  years 
include  similar  cases  where  downwelling  CTWs  trigger 
anticyclonic  eddies  east  and  southeast  of  the  GT.  These 
new  forcing  mechanisms  and  region  of  generation  (south¬ 
east  of  the  GT)  do  not  contradict,  but  complement  the  well- 


Figure  17.  SSH  anomaly  (color  contours  in  centimeters)  for  18  different  dates  in  June- August  1997  as  simulated  by 
1/16°  Pacific  NLOM  (1.5-layer  reduced  gravity  configuration).  The  positions  of  the  Gulf  of  Tehuantepec  (GT), 
Acapulco  (ACA),  Cabo  Corrientes  (CC),  and  Maria  Islands  (MI)  are  indicated. 
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Figure  18.  SSH  anomaly  (color  contours  in  centimeters)  for  1  8  different  dates  in  June- August  1997  as 
produced  by  the  Modular  Ocean  Data  Assimilation  System  [Fox  et  ai ,  2002].  The  positions  of  the  Gulf 
of  Tehuantepec  (GT),  Acapulco  (ACA),  Cabo  Corrientes  (CC),  and  Maria  Islands  (MI)  are  indicated. 
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accepted  theory  of  anticyclonic  Tehuantepec  eddies  gener¬ 
ated  by  strong  wind  events  in  the  western  GT. 

4.  Summary  and  Concluding  Remarks 

[29]  For  the  First  time  Tehuantepec  eddies  have  been 
studied  using  a  high-resolution  numerical  ocean  model 
forced  with  realistic  interannual  high-frequency  winds. 
Results  from  the  model  and  an  almost  8-year  time  series, 
October  1992  to  March  2000,  of  sea  surface  height  mea¬ 
sured  by  the  TOPEX/Poseidon  satellite  altimeter  show 
interannual  variability  in  the  number  and  strength  of  the 
Tehuantepec  eddies  (Figure  2  and  Table  1).  Demonstration 
of  the  Tehuantepec  eddy  interannual  variability  was  made 
possible  by  the  ability  of  the  satellite  altimeter  to  measure 
sea  surface  height  during  all  weather,  all  seasons,  and  all 
oceanic  conditions.  A  surprising  and  notable  finding  is  that 
the  altimeter  measurements  and  the  model  results  include 
anticyclonic  Tehuantepec  eddy  formation  during  a  5-month 
period  of  weak  Gulf  of  Tehuantepec  winds  throughout 
summer  of  1997  (Figures  2,  3,  13,  14,  17,  and  18),  and 
anticyclonic  eddy  generation  east  and  southeast  of  the  Gulf 
of  Tehuantepec  (Figures  5,  8,  13,  14,  17,  and  18).  Addi¬ 
tional  examples  of  eddy  generation  in  the  area  with  strong 
Gulf  of  Tehuantepec  winds  are  found  in  other  years  and 
other  seasons  of  the  same  year.  These  new  facts  comple¬ 
ment,  but  do  not  invalidate,  the  well-accepted  theory  of 
anticyclonic  Tehuantepec  eddy  formation  during  fall/winter 
on  the  western  side  of  the  Gulf  of  Tehuantepec  driven  by 
strong  Gulf  of  Tehuantepec  wind  events  (Figures  5  and  6). 
Equally  remarkable  is  that  the  model,  which  does  not 
assimilate  any  oceanic  data,  reproduces  qualitatively 
(Figure  2),  and  quantitatively  (correlation  coefficient  of 
0.72)  the  altimeter  sea  surface  height  observations. 

[30]  Opposite  to  our  initial  hypothesis,  Tehuantepec  eddy 
formation  and  Tehuantepec  eddy  interannual  variability 
cannot  be  explained  solely  in  terms  of  strong  and  intermit¬ 
tent  Tehuantepec  wind  events  (Figure  3).  In  fact,  two 
different  anticyclonic  Tehuantepec  eddy  generation  mecha¬ 
nisms  are  identified:  first,  the  generation  and  modulation  of 
the  Tehuantepec  eddies  solely  due  to  the  Tehuantepec  winds 
and,  second,  the  anticyclonic  Tehuantepec  eddies  generated 
by  free  propagating  baroclinic  intraseasonal  and  interannual 
downwelling  coastally  trapped  waves.  During  their  north¬ 
ward  propagation  and  while  propagating  along  the  non¬ 
meridional  coasts  of  Central  America  and  Mexico,  these 
waves  produce  a  strong  horizontal  and  vertical  shear  of  the 
horizontal  velocity  (Figure  11),  and  trigger  barotropic  and 
baroclinic  instabilities  in  the  flow  (Figures  12  and  15-16) 
generating  anticyclonic  Tehuantepec  eddies  (Figures  5,  13, 
14,  17,  and  18).  These  equatorial ly  generated  waves  have 
significant  interannual  variability  increasing  (decreasing) 
during  El  Nino  (La  Nina)  years  (Figures  7).  In  contrast, 
the  Tehuantepec  and  Papagayo  winds  do  not  show  a 
significant  interannual  variability  (Figure  3).  That  suggests 
a  dominant  role  for  baroclinic  downwelling  coastally 
trapped  waves  in  the  modulation  of  the  altimeter  measured 
and  model  simulated  Tehuantepec  eddy  interannual  vari¬ 
ability  (Figure  2). 
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